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ABSTRACT 

Future large-scale structure surveys of the Universe will aim to constrain the cosmo- 
logical model and the true nature of dark energy with unprecedented accuracy. In 
order for these surveys to achieve their designed goals, they will require predictions 
for the nonlinear matter power spectrum to sub-percent accuracy. Through the use of 
a large ensemble of cosmological ./V-body simulations, we demonstrate that if we do 
not understand the uncertainties associated with simulating structure formation, i.e. 
knowledge of the 'true' simulation parameters, and simply seek to marginalize over 
them, then the constraining power of such future surveys can be significantly reduced. 
However, for the parameters {n s , h, fif,, f2 m }, this effect can be largely mitigated by 
adding the information from a CMB experiment, like Planck. In contrast, for the 
amplitude of fluctuations as and the time-evolving equation of state of dark energy 
{wq, w a }, the mitigation is mild. On marginalizing over the simulation parameters, we 
find that the dark-energy figure of merit can be degraded by ~ 2. This is likely an 
optimistic assessment, since we do not take into account other important simulation 
parameters. A caveat is our assumption that the Hessian of the likelihood function 
does not vary significantly when moving from our adopted to the true simulation pa- 
rameter set. This paper therefore provides strong motivation for rigorous convergence 
testing of iV-body codes to meet the future challenges of precision cosmology. 

Key words: Cosmology: large-scale structure of Universe. 



1 INTRODUCTION 

Future spectro/imaging surveys of the low-redshift Universe 
such as DElfl KiDlfl EuclioQ and WFIRS10 will aim to 
constrain the cosmological model to unprecedented accu- 
racy. This will require impressive handling of every step of 
the observational pipeline in order to limit the possibility of 
systematic errors that may degrade the constraints on cos- 
mological parameters. Besides the observational processing, 
there will also be a similarly high demand placed on our abil- 
ity to generate theoretical predictions that are sufficiently 
accurate not to bias inferred cosmological parameters. These 
predictions will take the form of a set of estimators for the 
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primary observables that we intend to measure and their 
covariance, and also most likely their cross-covariance. The 
latter will be required for rob ustly testing for modifications 
to gravity (|Reves et alj|20ich . 

In galaxy clustering or cosmic shear surveys of the Uni- 
verse, the primary observables of interest are related to 
the matter power spectrum and its evolution with redshift. 
The matter power spectrum is the two-point covariance of 
the matter fluctuations transformed into Fourier space. The 
power spectrum provides a wealth of information on th e 
cosmological parameters (|Dod elson 2003; Wein berg! 120081 ). 
In order to maximize the amount of information obtainable 
from the power spectrum, we need to understand its de- 
pendence on the cosmological parameters in the nonlinear 
regime. 

A number of theoretical and semi-empirical techniques 
are available for predicting the nonlinea r power spectrum: 
such as the halo mode l dSeliakl l2000l : IPeacock fc Smith! 
l2000l : llvla fc Fry! |2000| : ISmith et all l2003l b perturba- 
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tion theory (|Bernardeau et aTll2002l ; 

2008) ; and scale tran sformations l|Hamilton et al,l 11991 



Crocce &: Scoccimarrol where the Fourier modes of the density field are given by 



Peacock fc Dod ds 1996). However, it is widely acknowledged 



that JV-body simulations provide the most direct path to- 
wards this answer. Cosmological A^-body s imula t ions are not 
witho ut pit-falls (|Heitmann et al.l 120051 . 120081 ; iReed et all 
2012), and in general they depend on a number of pseudo- 
free parameters, such as: the number of particles used to rep- 
resent the phase-space; the box-size; the redshift at which 
the initial conditions are given; the scale on which two-body 
forces are softened, etc. If a particle-mesh (PM) approach is 
employed then one additionally needs to set the scale above 
which forces are solved using mesh based techniques. If a 
tree technique is used, then one additionally needs to adopt 
a choice for the multipole order to which forces are expanded 
and the type and depth of the tree to be used. If both are 
used, then one also needs to set parameters that interpolate 
between the tree and PM methods. 

Given the complexity of the state-of-the-art TV-body 
codes, we are then lead to ask the following questions: 

• How do we determine the values of the true simulation 
parameters? 

• How much would forecasted parameter constraints be 
degraded if we were to marginalize over the simulation pa- 
rameters? 

• How does this affect the dark energy 'figure of merit'? 

In this paper we shall employ a large ensemble of N- 
body simulations to directly answer these latter two ques- 
tions, and leave the first for future study. 

The paper can be broken down as follows: In we 
provide a brief overview of the necessary theoretical back- 
ground, and in particular we give a brief review of the Fisher 
matrix approach to forecasting cosmological constraints. In 
[J3]we describe the large ensemble of simulations that we em- 
ploy. In i|4] we present results concerning the fiducial model 
power spectrum and its covariance matrix. In fJ5]we explore 
the dependence of the matter power spectrum on cosmo- 
logical and simulation parameters. In [|6]we use the Fisher 
matrix approach to explore how various assumptions con- 
cerning our understanding of iV-body simulations can im- 
pact the cosmological parameter forecasts from future large- 
scale-structure surveys. In Sj7| we focus on the question of 
constraining the time evolution of dark energy and how igno- 
rance of simulation parameters impacts the figure of merit. 
Finally, in SjH]we summarize our findings and conclude. 



2 FORECASTING COSMOLOGICAL 
CONSTRAINTS 

2.1 The Gemeinsam likelihood function 

We are interested in forecasting the ability of a future sur- 
vey of the universe to constrain the cosmological parameter 
space. We may assess this using the Fisher matrix approach. 

Consider a particular statistic that we will estimate 
from the survey data, and let us be concrete and take this to 
be the matter pow er spectrum P (k). The power spectrum 
may be defined as (|Peebleslll980l ): 



<5(k) = ^- f d 3 xexp [-ik ■ x] 5(x) 
and where the over-density field is defined as: 



5(x) 



P(x) - P 



(2) 



(3) 



with p and p being the local and background density. 

A given theoretical cosmological power spectrum de- 
pends on the wavenumber k - here we focus on the real-space 
isotropic function - and also the cosmological parameters 0. 
We are also interested in the case where the theoretical pre- 
dictions also depend on a set of internal simulation param- 
eters V- Let us write the augmented vector of cosmological 
and simulation parameters as <p — (6,ip). We denote the 
measurement of P(k) at wavenumber by Pi and the the- 
oretical (simulated) spectra by P slm (fc;|</>). Notice that here 
we are making the approximation that Pi does not depend 
on 8; this in fact is not true and any measurement of P re- 
quires the assumption of a cosmological model (we reserve 
further discussion of this for future work and note that this 
simply makes the analysis sub-optimal). 

Let us adopt a Bayesian approach to the analysis of our 
data and write the m measurements of the power spectra at 
wavenumbers k — » {ki, . . . , k m }, as P — > {Pi, . . . , P m }- The 
probability that our survey yields observations P, given the 
cosmological and simulation parameters <f>, is L(P|</>) - the 
likelihood. If the likelihood is Gaussian, then we have 



L(P|<£) 



• exp 



(4) 



(2tt)"/ 2 |C| 1 /2 

where we have made use of the Einstein summation conven- 
tion. In the above equation we also defined 

-yrthcory 



Vi(4>) 



'fa]*) 



(5) 



where P the ° ry (fci|</>) is the expectation for the theory power 
spectrum. C(</>) is the covariance matrix, which may be de- 
fined as: 



p theory (fc, 



(6) 



and |C| is the determinant of the covariance matrix. 

Using Bayes theorem, the likelihood is directly related 
to the posterior probability, p(<p\P), through a set of priors, 
p{4>), an d is normalized by the evidence, p(P): 



pWP) = 



p(cp)L(P\ 
p(P) 



p(<t>)L(P\4>) 

f d<PP(<p)L(P\<t>) 



(7) 



If the priors are flat, then the posterior probability is simply 
proportional to the likelihood. Close to its maximum, at <j>o, 
we may Taylor expand the logarithm of the posterior, and 
for flat priors also the log likelihood {C = lnL), to obtain: 



Inp(0[P) 



£(P|0) 
£(P|0o) 



-&.(j) a 'H a fj{<pQ)&-4>p + 



(8) 



where in the above: A<j> = (</> — 0o) are deviations of the 
parameters from the fiducial values; the first derivative van- 
ished at the maximum; the second derivative is identified 

as 



U M <<5(ki)<5(k 2 )> = P(fci)C- 



k 2 



(1) 



T-La0 



(9) 
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and it is given the name of Hessian, or curvature matrix. 
We have also used the notation: = d/d<j) a - In truncating 
this expression for the posterior at second order we are im- 
plicitly assuming that the likelihood is also Gaussian in the 
parameters. Hence, we may rewrite the above expression for 
the posterior as, 



K0) 
P (P) 



L(4>o) exp 



'-A(f) a 'Hap{4 , o)^4'0 



(10) 



Thus Hep informs us about errors on the parameters and 
how different parameters may be correlated with respect to 
each other - in the context of their effects on the data. 

Since the likelihood itself depends on the data, it is also 
a random variable. Taking an ensemble average over many 
realizations of the data, we arrive at the Fisher matrix: 



Tap = (7~L a p) = — 



d 2 lnL \ 

d(f> a dct>i3 / 



(11) 



Considering the division into cosmological and simula- 
tion parameters, this matrix may be written schematically 
as: 



r 



■4><t> 



(12) 



where T 66 ', 7^ and T i>M> denote the Fisher matrices of the 
cosmological, cosmological-cross-simulation and simulation 
parameter spaces, respectively. 

From the Fisher matrix one can obtain the expected 
marginalized error on parameter <j>i and the covariance be- 
tween parameters (<j)i, cj>j): 



4"l>] 



<W1 



(13) 



We can also obtain conditional errors for the cosmological 
parameters, conditioned on the simulation parameters pos- 
sessing a particular value: 



[J* 



(14) 



These expressions represent the minimum variance bounds 
(MVB) (for a derivation see lHeavenj|2009l ). 

Lastly, we note that for the specific case of a Gaussian 
likelihood, it can be shown that the F i sher matrix ta kes on 
the special form dTegmark et al.lll997l ; lHeavensll2009l '): 



F a0 = -Tr [C- 1 C,„C- 1 C,p] 



^theory 



c -l p thco ry ^ (15) 



Our first objective may now be reformulated as the fol- 
lowing questions: 

• How do the MVBs obtained from T 9B compare with 
those for T^^l How do our parameter forecasts degrade 
when we marginalize over simulation parameters? 



3 THE iV-BODY SIMULATIONS 

In order to study the Fisher information we have generated a 
large suite of iV-body simulations. As Eq. (|15[) demonstrates, 
one needs to compute the derivatives of the theoretical power 
spectra with respect to the parameters cj> and also the co- 
variance matrix. In fact, one also needs the derivatives of the 
covariance matrix with respect to the parameters (p. Since 



Table 1. Cosmological parameters used for the fiducial zHORIZON 
suite of simulations and the variations with respect to the cos- 
mological parameters. The columns are: name of simulation se- 
ries; density parameters for matter, dark energy and baryons; the 
equation of state parameter for the dark energy P w = uip w ; nor- 
malization and primordial spectral index of the power spectrum; 
dimensionless Hubble parameter h, respectively. 



Par am. 




Q,, 


U'() 




(Xg 


n 


h 


Fid. 


0.25 


0.04 


-1.0 


0.0 


0.8 


1.00 


0.7 


VI 


0.20 


0.04 


-1.0 


0.0 


0.8 


1.00 


0.7 


V2 


0.30 


0.04 


-1.0 


0.0 


0.8 


1.00 


0.7 


V3 


0.25 


0.035 


-1.0 


0.0 


0.8 


1.00 


0.7 


V4 


0.25 


0.045 


-1.0 


0.0 


0.8 


1.00 


0.7 


V5 


0.25 


0.04 


-1.2 


0.0 


0.8 


1.00 


0.7 


V6 


0.25 


0.04 


-0.8 


0.0 


0.8 


1.00 


0.7 


V7 


0.25 


0.04 


-1.0 


-0.1 


0.8 


1.00 


0.7 


V8 


0.25 


0.04 


-1.0 


0.1 


0.8 


1.00 


0.7 


V9 


0.25 


0.04 


-1.0 


0.0 


0.7 


1.00 


0.7 


V10 


0.25 


0.04 


-1.0 


0.0 


0.9 


1.00 


0.7 


Vll 


0.25 


0.04 


-1.0 


0.0 


0.8 


0.95 


0.7 


V12 


0.25 


0.04 


-1.0 


0.0 


0.8 


1.05 


0.7 


V13 


0.25 


0.04 


-1.0 


0.0 


0.8 


1.00 


0.65 


V14 


0.25 


0.04 


-1.0 


0.0 


0.8 


1.00 


0.75 



Table 2. The Gadget-2 parameters used for all fiducial simu- 
lations, and the variations with respect to the simulation pa- 
rameters. The columns are: name of simulation series; simula- 
tion parameter varied; fiducial value; simulated variations. Pa- 
rameters that were not varied include: MinSizeTimestep=0.000; 
TypeOf 0peningCriterion=l; TreeDomainUpdateFrequency=0.05; 
PartAllocFactor=2.0; TreeAllocFactor=1.0 



Simulation 


Parameters 


Fiducial 


Low 


High 


S1/S2 


ErrTolForceAcc 


0.005 


0.004 


0.006 


S3/S4 


ErrTolIntAcc 


0.025 


0.02 


0.03 


S5/S6 


ErrTolTheta 


0.45 


0.4 


0.5 


S7/S8 


PMGRID 


750 


500 


1000 


S9/S10 


MaxRMSDispFac 


0.2 


0.15 


0.25 


S11/S12 


Softening 


0.03 


0.025 


0.035 


S13/S14 


RCUT 


4.5 


4.0 


5.0 


S15/S16 


ASMTH 


1.25 


1.15 


1.35 


S17/S18 


MaxSizeTstep 


0.025 


0.020 


0.03 



estimating the covariance matrix is a sufficiently challenging 
task in itself, we shall reserve the inclusion of information 
from this for future study. Henceforth, we shall drop the first 
term in Eq. (|15[) from ou r analysis (for f urther justification 
of this approximation see lTegmarkHl997T ). 

In order to determine the covariance matrix we 
have simulated 200 realizations of our fiducial cosmolog- 
ical model. The specific cosmological parameters that we 
adopted are for a flat ACDM model with: {as = 0.8, fi m = 
0.25, fi 6 = 0.04,^o = -1.0, w a = 0.0, h = 0.7, n s = 1.0} 
where: erg is the variance of mass fluctuations in a top-hat 
sphere of radius R = 8/i _1 Mpc; fl m and fij, are the matter 
and baryon density parameters; wo and w a are the constant 
and time-evolving equation-of-state parameters for the dark 
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energy, i.e. Pde/pde = w(a) — wo + (1 — a)w a ; h is the di- 
mensionless Hubble parameter; and n is the power-law index 
of the primordial density power spectrum. Our adopted val- 
ues were inspired by the results of the WMAP experiment 
(|Komatsu et alj|2009l ). Table [T] contains further details of 
the cosmological parameters of the fiducial model. 

All of the iV-body simulations were run on the zBOX-3 
and Schrodinger supercomputers at the University of 
Zurich, using th e publicly available Tree-PM code GADGET-2 
(|Springeli r2005). with a slight modification that permitted 
a time-evolving equation of state for dark energy, specified 
by the parameters {wo,w a }- This code was used to follow 
with high-force accuracy the nonlinear evolution under grav- 
ity of N = 750 3 equal mass particles in a comoving cube 
of length L = 1500 /i _1 Mpc, giving a total sample volume 
of order V ~ 540 h~ 3 Gpc 3 . Newtonian two-body forces are 
softened below scales i S oft = 60 /i _1 kpc. We shall refer to 
this suite of simulations as the zHORIZON runs (Zurich Hori- 
zon simulations). The transfer function for the simulations 
was generated using the p ublicly available cmbfast code 
(|Seliak fc Zaldarriaga|[l99a ). with high sampling of the spa- 
tial frequencies on large scales. For the time evolvin g dark 
energy models we used the code C AMB ([Lewis et al.i 2000) 
and with the dark energy module of iHu fe Sawickil (|2007f ). 



Initial conditions were lain down at redshift z = 49 
using the serial ver sion of the publicly available 2LPT code 
(|Crocce et alj|2006l ). The zBOX-3 runs took roughly ~20Hrs 
per run on 256 cores, and the Schrodinger runs took ~6Hrs 
per run on 256 cores. For all of the realizations snapshots 
were output at a number of redshifts, though for this study 
we focus only on the results at z = {1, 0.5, 0}. For complete- 
ness, the Gadget-2 parameters that we used are presented 
in Table H 

In order to evaluate the derivatives of the power spec- 
trum with respect to the cosmological parameters, we have 
performed an additional 56 simulations - the cosmological 
variations, labeled V1-V14. We have considered the effect 
of changing a single cosmological parameter, whilst holding 
the remaining parameters fixed. For each such modification 
we ran 4 simulations. We used double-sided variations, e.g. 
P(k\4> + A(f>i) and P(k\(j> — Acj>i), as this enables more accu- 
rate computations of the numerical derivatives, which will 
be important for our Fisher-matrix estimates. Also, in or- 
der to reduce the noise in these estimates, we matched the 
initial Gaussian random field of each realization with the 
corresponding one from the fiducial model. Full details of 
these simulations are summarized in Table [T] 

To estimate the derivatives of the power spectrum with 
respect to the simulation parameters, we have performed 
another 18 simulations - the simulation variations, labeled 
S1-S18. This time we keep the cosmological parameters as 
the fiducial ones and explore the effect of changing a sin- 
gle simulation parameter, whilst holding the remaining ones 
fixed. For each such modification we ran a single simula- 
tion, but again we considered double-sided variations, with 
matched initial Gaussian random fields so as to decrease 
the noise when estimating derivatives. The exact list of sim- 
ulation parameters that we have sampled are presented in 
Tabled 




k [h/Mpc] 

Figure 1. The ensemble- averaged dark matter power spectrum 
for the 200 fiducial realizations, at redshifts z E {1,0.5,0}, with 
errors per realization. Top panel: absolute dimensionless power, 
A 2 = k i P(k)/2-K 2 . Central panel: the ratio of the mean power 
spectra with respect to the linear theory. Bottom panel: ratio 
of power spectra with respe ct to the nonlinea r predictions from 
the fitting formula halof it l lSmith et al,ll2003l ). In all panels the 
solid and dashed lines denote the nonlinear and linear theory 
predictions, respectively. Note that in the central and bottom 
panels, the z = 0.5 and z = results have been off-set by 0.1 and 
0.2 in the vertical direction for clarity. 



4 ANALYSIS I: THE FIDUCIAL MODEL 
4.1 Power spectrum 

As a first exploration of the simulation data we compute the 
matter power spectrum at the redshifts of interest. 

The power spectrum in the simulation cube for a given 
Fourier mode is as described in Eq. JT}. In practice, the 
power is estimated by averaging over all wavemodes in thin 
spherical shells in fc-space - band-powers. The band-power- 
averaged power spectrum can be written, 

U 1 d a k (s i Qc)S i (-k) S 



Vs,i J V 

= ^E( 5<i ( k ^ d (-M) , (16) 

* 3=1 

where the average is over the k-space shell V 3 , of volume 

r ki + Ak/2 



V s ,i = 



d 3 k = AirkfAk 



-Afc/2 



+ 12 V h , 



(17) 



© 0000 RAS, MNRAS 000, 000-000 



Cosmological constraints and N-body codes 5 



and where = V s ,i/Vk is the total number of modes in 
the shell. Vk = kf is the fundamental fc-space cell volume 
and kf = 2n/L is the fundamental wavemode. 

Notice that in Eq. (|16|l we have used the superscript d, 
this stands for discrete, since we make a Fourier decompo- 
sition of the point-sampled field. The point-sampled power 
spectrum is related to tha t of the continuous mass d ensity 
field through the relation l|Peebleslll980l ; ISmithll2009l ) : 

P d (k) = P c (k) + - , (18) 
n 

where P c is the power spectrum of the underlying continu- 
ous field. The constant term on the right-hand side of the 
equation is more commonly referred to as the 'shot-noise 
correction' term, and is the additional variance introduced 
through discreteness, where n = N/Vfj,. 

In order to compute the power spectrum, we apply the 
standard F ast Fourier Transform (FFT)-based approach (se e 
for example [Smith et al.ll20oi|jindl2005l ; ISmith et al.ll2008l l. 
We use a 'cloud- in- cell' (CIC) mass assignment scheme for 
the simulation particles, and deconvolve each Fourier mode 
accordingly with the window function. The power spectrum 
estimator is then given by Eq. (|16[) . We use FFT grids with 
A'grid = 1024 cells per dimension, and this sets the mini- 
mum and maximum spatial frequencies to: k m i n = 2n/L — 
0.0042 /iMpc" 1 and fc Ny = TiN gM /L = 2.15 ftMpc" 1 . In 
practice, the power on length scales k > &ny will get 
'aliased' to larger spatial scales, and so we take fc ma x = 
fcNy/2. In this study we have decided to estimate the power 
spectrum in 35 logarithmically spaced band powers in the 
interval k G [0.0042, 1.0] ftMpc -1 . We adopt this strategy in 
order to obtain sufficiently high signal-to-noise estimates of 
the covariance matrix. 

Figure [T] shows the ensemble-averaged dark matter 
power spectrum for the 200 realizations, at the redshifts 
z G {1,0.5,0}, denoted by the green, red and blue points, 
respectively. The coloured points actually show the power 
spectra obtained from a linear binning, where the bin spac- 
ing is in units of the fundamental fc-cell spacing, kf. The 
black points denote the results for the 35 logarithmically 
spaced bins, and the error bars are on the mean. In the top 
panel of the figure we show the dimensionless power, which 
may be defined: 

A 2 (fc) = J^s^P(k) . (19) 

In the middle panel we show the ratio of the measured power 
spectra with respect to the input linear theory power spec- 
tra. For clarity, we offset the power spectra at z = by 
20% and at z = 0.5 by 10% in the vertical direction. The 
black solid line show s the nonlinear po wer spectra predic- 
tions from halofit (|Smith et al. 1 120031 '). The plot demon- 
strates that strong nonlinear amplification occurs on scales 
k > 0.1 ftMpc -1 , and that linear theory is not a good ap- 
proximation on these scales. In addition, nonlinear ampli- 
fication is not significantly weaker at higher redshifts. The 
bottom panel shows the ratio of the measured power spec- 
tra with respect to halofit. Again we have offset the power 
spectra for clarity, halofit is able to describe the measured 
power spectra to better than 10% on the scales investigated. 

The shot-noise correction to the power spectrum is 
fshot = 1/n = 8ft _3 Mpc 3 , which in terms of the dimen- 
sionless power is Ag hot = 4fc 3 /7r 2 . At k ~ 1/iMpc -1 , this is 




0.01 0.1 1 

k [h/Mpc] 

Figure 2. Ratio of the measured error in the power spectrum 
to the Gaussian-predicted error, i.e. Eq. 1241 The measured error 
is obtained from the ensemble of 200 JV-body simulations. Blue, 
green and red points depict the results for redshifts z G {0, 0.5, 1}, 
respectively. The solid and open symbols show the results ob- 
tained when the linear theory and measured nonlinear power 
spectra are used to make the predictions from the Gaussian error 
formula. For clarity, the z = and z = 0.5 results have been 
off-set by factors of 3 and 2 on the y-axis, respectively. 

A 2 hot ~ 0.4. Thus for z = 1, the power has roughly a 10% 
correction, which by z = is reduced to 2%. 

4.2 Covariance matrix 

An unbiased estimator for the covariance between different 
band power estimates can be obtained through: 

1 n e 

e « = jv^tEP"'-^!?'-^ ; (2°) 

CH — 1 

* = (2D 

a — 1 

where JVe is t he number of realiz a tions. 

Following IScoccimarro et al.1 (ll999T ) and ISmitbl ||2009h . 
a theoretical expression for the bin-averaged covariance ma- 
trix of the matter power spectrum, obtained from a set of 
densely-sampled tracers of the mass field, can be written: 

c -=^ + wM 6 ^' (22) 

where Ty is the shell-averaged, connected part of the trispec- 
trum in parallelogram configuration: 

Ti, = [ f. kld ' k2 f(k 1 ,k 2 ,-k 1 ,-k 2 ) , (23) 

J Vs,i *s,j 

with f(k 1 ,k 2 ,k 3 ,k 4 )V^ = <<5(k 1 )...5(k 4 )}<5£ + ... +k4 , be- 
ing the matter trispectrum. Note that for a Gaussian ran- 
dom field the connected part of the trispectrum vanishes, 
i.e. T = 0, and the covariance reduces to: 



© 0000 RAS, MNRAS 000, 000-000 



6 R. E. Smith et al. 




0.01 0.1 1 0.01 0.1 1 0.01 0.1 1 

k [h/Mpc] 



Figure 3. Top panel: Evolution of the power spectrum correlation matrix as a function of wavenumber, estimated from the ensemble of 
200 simulations. The left, central and right panels show the results for z = {1.0, 0.5, 0.0}, respectively. Bottom panel: same as above, 
only the correlation matrix has been box-car smoothed by a square top-hat filter of size 3x3 pixels. This clearly reduces the noise in 
the correlation matrix on large scales. 



Figure [5] shows the standard deviation for the matter 

1/2 

power spectrum, i.e. , measured from the simulations, 
scaled in units of the square root of the Gaussian expecta- 
tion for the variance given in the equation above. The figure 
reveals that for k < 0.1/iMpc" 1 the diagonal errors are 
reasonably well described by Eq. (|24[) . However, on smaller 
scales we find that the errors are significantly larger than 
one would expect from simple mode counting. If one uses 
Eq. (|24[) with P measured from the simulations (open points 
in Fig.0), then the errors appear to increase oc k. This would 
suggest that in the deeply nonlinear regime T« oc Pf as 
k — > oo. This scaling is consisten t with the predictions fro m 
the 1-Loop perturbation theory l|Scoccimarro et alJll999h . 

An interesting way to present the information in the 
off-diagonal elements of the covariance matrix is through 
the cross-correlation matrix. This may be defined as: 

r« = Cij , (25) 

V CiiCjj 

and it gives the strength of the covariance in a particular el- 
ement relative to the square root of the product of variances 
in the relevant bins. It is therefore bound to the interval 
re [-1,1]. 

The top panels in Figure [3] present the evolution of the 
correlation matrix as a function of redshift from z = 1 to 
2 = 0. The increasingly redder/bluer colours demonstrate 
increasing/decreasing correlation strength. The results on 
large scales appear to be slightly noisy. The bottom panels 
of Fig. [3] presents the same information, only here we have 



performed a box-car smoothing of the correlation matrix in 
order to reduce the noise. For each pixel, we take the av- 
erage of all pixels that are within 1-pixel from the current 
pixel centre, excluding the pixels on the diagonal and being 
careful in our treatment of the edges of the matrix. We also 
keep the diagonal elements fixed at r = 1 (for furth er dis- 
cussion of this approach see iMandelbaum et al.ll2012h . This 
noise reduction strategy constitutes a plausible alternative 
to various oth er ad-hoc approa ches presented elsewhere in 
the literature (|Ngan et alJl2012T l. 

For the case of both the raw and the noise-reduced ma- 
trix, the off-diagonal correlations are in general non-zero and 
positive. The correlations increase as the wavenumbers of 
the two considered band powers are increased. Also, the cor- 
relation increases with decreasing redshift. For our choice of 
binning and simulation volume, we find that different power 
spectral band powers are > 50% correlated for {ki,kj} > 
0.2/iMpc" 1 at z = 0, and for {ki,k 3 } > 0.25 hMpc' 1 by 
z = l. 

Figure [4] presents slices through the power spectrum 
correlation matrices measured at z £ {0,0.5, 1} for the 200 
realizations of the fiducial model. These results point to a 
reasonably good agreement between the raw and box-car- 
smoothed covariance matrices. 

The covariance matrix of the matter pow e r spec trum 
has recently been explored by iTakahashi et alj (|2009l ) who 
ran 5000 PM simulations in boxes of size L = l/i _1 Gpc 
with N = 256 3 particles. These simulations are not of 
sufficiently-high spatial resolution to probe the covariance 
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Figure 4. Different rows of the power spectrum correlation matrix as a function of wavenumber. From left to right the columns show 
results for k\ = {0.01, 0.02, 0.04, 0.12, 0.32, 0.84} h Mpc -1 , as a function of ki- The top, middle and bottom rows show the results for 
z = {0, 0.5, 1}, respectively. The dashed and solid lines show the estimates for the raw and box-car smoothed correlation matrix. 



of power spectrum estimates beyond scales of the order 
k ~ 0.2(0.4) /iMpc" 1 at ~ 1(3)% precision. Moreover, whilst 
they did employ the more accurate 2LPT initial conditions 
- as does our study - they also used a rather low start red- 
shift of z = 20, w hich may induce small scale inaccuracies 
(|Reed et al.ll2012h . On comparing their results with ours, we 
note that whilst they have a factor of 25 times more simu- 
lations, each of our simulations has a factor of 3 times more 
volume. This makes the overall difference roughly a factor 
of ~ 3 in terms of (S/N). We should therefore be able to ob- 
tain a reasonably accurate covariance matrix. This is further 
mitigated by our smoothing of the correlation matrix. 

W e also compare our study with that of iNgan et alj 
<|2012h . who used the code CUBEP3M to run 1000 simula- 
tions of boxes with L = 600 /i" 1 Mpc and with N = 256 3 
particles. We underline that for this choice of simulation 
set-up, the shot-noise corrections to the power spectrum at 
k — 1/tMpc -1 are 6% and 30% at z = and z = 1, re- 
spectively. Again, whilst their study used 1000 simulations, 
our simulations have roughly 15 times more volume per run. 
Moreover, they have explored the covariance matrix for 54 
bins, nearly a factor of 2 times more than we employ, hence 
the relative statistical power of our study should be at the 
very least co mparable w i th the ir work. It is also worth point- 
ing out that lNgan et alj (l2012h found a 20% anti-correlation 
of band-powers on the largest scales. We find no evidence of 
such a strong anti-correlation. We also note t hat th e power 
spectra from the simulations of INgan et al l (|2012h appear 
to show a worrying ~ 5 — 7% positive off-set from the linear 
theory predictions for the power spectrum on scales com- 
parable t o the box-scale, although t his issue may now be 
resolved l|Harnois-Deraps et al.ll2012h . 

A more recent study bv lde Putter et all (|2012h has used 



a suite of 1024 simulations of L = 600 h 1 Mpc boxes and 
160 simulations of a L = 2400 /i _1 Mpc box to explore the co- 
variance matrix. They found that the results from the large- 
box simulations were in reasonably good agreement with the 
larger ensemble of smaller box simulations. 

Thus our results are in broad agreement with all of these 
works and non-trivial band power correlations must be ac- 
counted for in cosmological analysis of large-scale structure 
data. 



5 ANALYSIS II: VARIATIONS 

5.1 Power spectrum dependence on cosmological 
parameters 

We now turn to the task of exploring the cosmological de- 
pendence of the power spectrum. As mentioned in fj3] we 
consider the variations with respect to 7 cosmological pa- 
rameters: = {(78, O ra , Qi, Wo, w a , h, n s }. For each cosmo- 
logical parameter, we freeze all of the other parameters and 
simulate two variations, up and down, around the fiducial- 
model value. For each such variation we have performed 4 
realizations. 

Figures [BT ll 1 1 present the ensemble-averaged variations. 
For each figure the left, middle and right panels show the 
results at redshifts: z = {1, 0.5, 0}, respectively. The top sec- 
tions show the absolute power; the middle, the ratio with re- 
spect to the linear theory of the fiducial model; and the bot- 
tom the ratio with respect to the measured nonlinear power 
in the fiducial model. In all panels the red points denote the 
fiducial model, the blue stars denote the upper variation, i.e. 
P(ki\4> + A<ji M ), and the open green stars denote the lower 
variation, i.e. P(ki\tj> — A(/>^). It is worth noting that when 
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Figure 5. Dependence of the nonlinear matter power spectrum on the power spectrum normalization parameter erg. In all panels the 
solid blue and open green stars depict the estimates from the ./V-body simulations with as = {0.9, 0.7}, and the solid red points denote 
the results for the fiducial model <jg = 0.8. The top panels show the absolute power spectrum; the central panels show the ratio of the 
spectra with respect to the fiducial linear theory predictions; the bottom panels show the ratio of the spectra with respect to the fiducial 
power spectrum. From left to right, the three columns represent results for epochs z = {0, 0.5, 1}, respectively. 
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Figure 6. The same as Fig. [5] only this time showing the dependence of the power spectrum on the matter density parameter O m . The 
solid blue and open green stars depict the results for fl m = {0.3, 0.2}, and the solid red points denote the results for the fiducial model 
fi m = 0.25. 
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Figure 7. The same as Fig. [5] only this time showing the dependence of the power spectrum on the baryon density parameter Q.^. The 
solid blue and open green stars depict the results for Qf, = {0.045,0.035}, and the solid red points denote the results for the fiducial 
model Q m = 0.04. 
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Figure 8. The same as Fig. [5] this time showing the dependence of the power spectrum on the dark energy equation-of-state parameter 
wo- The solid blue and open green stars depict the results for u>o = {—0.8,-1.2}, and the solid red points denote the results for the 
fiducial model wq = —1.0. 
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Figure 9. The same as Fig. [5] this time showing the dependence of the power spectrum on the dimensionless Hubble parameter h. The 
solid blue and open green stars depict the results for h = {0.75, 0.65}, and the solid red points denote the results for the fiducial model 
h = 0.7. 
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Figure 10. The same as Fig. [5] this time showing the dependence of the power spectrum on the primordial-power-spectrum index n s . 
The solid blue and open green stars depict the results for n s = {1.05,0.95}, and the solid red points denote the results for the fiducial 
model n B = 1.0. 
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Figure 11. The same as Fig.[5] this time showing the dependence of the power spectrum on the dark energy equation-of-state parameter 
Wa- The solid blue and open green stars depict the results for w a = {0.1, —0.1}, and the solid red points denote the results for the fiducial 
model w a = 0. 



we compute the ratio of the variant power spectra with the 
nonlinear fiducial spectrum, we compute this ratio for each 
realization, before averaging. This leads to the cancellation 
of some of the cosmic variance, and explains why the error 
bars in the lower sections of each panel are not visible. 

We next turn our attention to the computation of the 
derivatives of the power spectra with respect to the cosmo- 
logical parameters. In order to obtain low-noise estimates of 
these, we take advantage of the matched initial conditions 
between the variations and use the double-sided derivative 
estimator, 



dP{k\4>) 



P{k\4>)- 



>logP(fe|<£) 



(26) 



where for the first term on the right-hand-side we take all 
200 of the fiducial simulations as described by Eq. (|16p . For 
the logarithmic derivative we use the estimator: 



aio g p(fc I |</») 



jV var 

cnscmb 



P (a) {h\4> + A«?v) - P (a) (h\4> - A<fv) 



2A0 /i P(°O(fc i |0) 



(27) 



cnsemb 

Figure [12] presents the evolution of the logarithmic 
derivatives of the power spectrum with respect to the 7 cos- 
mological parameters that we have considered. The deriva- 
tives are computed as described by Eq. (|27|l . In each panel 
the linear-theory derivatives are given by the solid blue lines 
and the black dashed lin es show the pred ictions for the 
derivatives using halofit (|Smith et al.ll200ot ). 



The figure demonstrates that on scales k < 0.1 feMpc -1 
one may capture the cosmological parameter dependence of 
the matter power spectrum through the variations in the 
linear power spectrum. However, on smaller scales on must 
obtain the derivatives from full nonlinear modelling. With 
the exception of wo and w a at z = 0, the predictions from 
halofit are in reasonably good agreement with the esti- 
mates from the numerical simulations. For wo and w a at 
z — halofit fails to predict the nonlinear derivatives. 
This may partially be explained by the fact that we have 
normalized the initial power spectra to have the same as'- 
had we instead adopted A„, the amplitude of the primordial 
power spectrum, as our power spectrum normalization cri- 
terion, then we expect that halofit would h ave made more 
reasonable predictions I Jennings et al.ll2010h . 

We also point out that as k — ¥ 1/iMpc -1 , the mea- 
sured derivatives for {fi m , fib, n s , h} appear to approach 
d log P/da — > 0. This suggests that there is very little 
cosmological information to be gained by the inclusion of 
measurements on very small scales. On the other hand, in- 
cluding the information from small scales can greatly in- 
crease the cosmological information about the parameters 
{as, wo, w a }- 



5.2 Power spectrum dependence on simulation 
parameters 

We now explore the dependence of the matter power spec- 
trum on the Gadget-2 simulation parameters. As described 
in fj3] we have considered variations in 9 of the parameters. 
As for the variations in the cosmological parameters, we take 
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Figure 12. Evolution of the logarithmic derivatives of the power spectrum with respect to the cosmological parameters. In all panels: 
solid red points denote the estimates from the ./V-body simulations; solid blue and black dashed lines denote the predictions from the linear 
theory and nonlinear halof it fitting function. From top to bottom the different rows show the results for variations in the parameters: 
{Q m , ag, n s , wq, h, Of,}; respectively. The left, central and right columns show the results for epochs, z = {0,0.5, 1}, respectively. 
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Figure 13. 

[P(fc|'0fid ± Aip) ~ P(k\ipfid)] /P{k\4>&d)- Blue, green and red points are for z 
points are for the positive and negative changes in the fiducial parameters. 
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a function of wavenumber. We plot 
{0, 0.5, 1}, respectively. The open and solid 



upper and lower variations of a single simulation parameter 
and freeze all others at their fiducial values. 

Figure Q3] presents the percentage differences between 
the variational and the fiducial models. Each of the nine 
panels corresponds to one parameter, with the solid and 
open points depicting the upper and lower variations respec- 
tively. The blue, green and red coloured symbols show the 
results at epochs z = {f,0.5, 0}, respectively. We find that 
the most significant source of error is given by the param- 
eter PMGRID, which can yield percent-level errors on small 
scales. We also find that the parameters which control the 
interpolation between the Tree- and the PM-force calcula- 
tions, RCUT and ASMTH, can also introduce significant, but 
sub-percent errors. Again, these are most important on small 
scales. ErrTolForceAcc and the Softening can also induce 
~ 0.2% errors in the power spectrum. The figure also shows 
that these variations are also relatively time-independent. 
This can be demonstrated more clearly by considering the 
logarithmic derivatives. 

Figure [T4] presents the logarithmic derivatives of the 
matter power spectrum with respect to variations in the 
Gadget-2 simulation parameters. We estimate the deriva- 
tives as described by Eq. (|27p , except that we only use a 
single realization to do this. We now make some important 
observations: firstly, on large scales, with the exception of 
the parameter ErrTolIntAcc, all of the derivatives are very 



close to zero. Moreover, they display a very weak depen- 
dence on redshift, which is a marked difference from the 
cosmological parameters, which tend to evolve with both 
time and scale. This is an important point, suggesting that 
the information coming from the cosmology dependence of 
the simulations can be disentangled from that coming from 
the simulation parameters. 



6 FISHER MATRIX RESULTS 

Having obtained estimates for the time evolution of the fidu- 
cial model power spectrum, its covariance matrix and its 
derivatives with respect to the cosmological and simulation 
parameters, we are now in a position to explore the true 
cosmological information content of the power spectrum. 



6.1 Estimator for the Fisher matrix 

We compute the Fisher matrix as described by Eq. (|f 5|l , and 
after dropping the first term, the estimator is: 

~ x - gjogg pr^i p d\ogP 3 , , 

- 2^ ^-fteT ■ (28) 
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Figure 14. Logarithmic derivatives of the power spectra with respect to the simulation parameters as a function of wavenumber. Blue, 
green and red points are for z = {0, 0.5, 1}, respectively. 



We shall not use the above equation directly, but an alter- 
nate form. Consider the matrix C, and let us rewrite it as 



Cij — CTi CTj Tij 



(29) 



(no summing over repeated indices), where erf is the variance 
associated to the ith measurement bin. The inverse of C can 
be written as 

C'j 1 = r^/ener-j . (30) 
Using the above identity allows us to rewrite Eq. (|28|l as: 



d log y t 

d<f>o 



~Vi r ijVj" 



dloe 



(31) 



where % = y{ki\ej>) = P(k i \(/))/er(k i ) and a(fej) = C'j/ 2 . This 
latter form is very useful, since one simply needs to invert 
the correlation matrix rather than the covariance matrix. 
In theory there should be no difference between the results 
from these two approaches, however on a computer there can 
be. Whilst the elements of the covariance matrix can differ 
wildly, even by orders of magnitude, the elements of the cor- 
relation matrix are constrained to range from [—1, 1]. Thus 
we reduce the risk of inaccurate and potentially unstable 
inverse estimates due to round-off errors. This is especially 
true when large dynamic ranges are considered and when 
many matrix elements are employed. 



6.2 Information content of the power spectrum 

For our cosmological forecast we adopt a survey consisting 
of three independent volumes, each of which has the same 
volume V M = 3.25 h~ s Gpc 3 , but mapping the three epochs 
z — {0, 0.5, 1}. Whilst this does not directly match a par- 
ticular survey, it covers the scale and evolution that should 
be obtainable with BOSS or DES. 

Figure [15] shows the fractional errors in the matter 
power spectrum for the variations in the 7 cosmological pa- 
rameters, at epochs z = {0, 0.5, 1}, and as a function of the 
maximum wavenumber that enters the calculation, i.e. fc max . 
Note that since the fiducial value of w a = 0, for this case we 
simply plot Aw a - In the figure, the unmarginalized errors 
on the parameters are given by the dotted black lines, i.e. 
Ap Q = J- aa - This 1-ct error is valid only if all the other 
parameters are known. 

If we are required to estimate all parameters from the 
data then the best we can ever do is to saturate the MVB, as 
described in !f2] If we assume that the simulation parameters 
are known, then the errors are given by Eq. (|13[1. i.e. Ap a = 



V \^ ee \aoii ant ^ we obtain the solid red lines in Fig. 1151 Notice 
that the cosmological information is significantly reduced. 
Once k ~ 0.4/iMpc -1 is reached, adding smaller scales does 
not reduce the errors on most of the parameters. This is with 
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Figure 15. Forecasted fractional \-a errors on the cosmological parameters as a function of the maximum wavenumber considered. In 
all panels, the dashed black lines denote the unmarginalized errors; solid red lines denote the errors marginalized over all remaining 
cosmological parameters; the blue dashed lines denote the errors after marginalizing over all other cosmological parameters and all 
simulation parameters. The top left through to top right panels present the results for a = {as , wo, Qh} and the bottom left to bottom 
right panels present the results for a = {O m ,n s , h,w a }. Note that for w a we simply plot Aw a . 
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Figure 16. Forecasted likelihood contours for the 15 possible parameter combinations. In all of the panels the solid red ellipses denote 
the \-a likelihood surface one expects from our fiducial survey marginalized over all other cosmological parameters. The blue dashed 
ellipses denote the same, but this time marginalizing over all other cosmological and simulation parameters. 
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Figure 17. The same as Fig. 1151 only this time we have added the information from a CMB experiment like the Planck satellite. Note 
that we have also employed a strong prior on the flatness of the Universe. 
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Figure 18. The same as Fig. 1161 only this time we have added the information from a CMB experiment like the Planck satellite. Note 
that we have also employed a strong prior on the flatness of the Universe. 
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the exception of wq, w a and as, for which adding small-scale 
structures does help. 

On the other hand, if we are to marginalize over the sim- 
ulation parameters, owing to the fact that we are ignorant 
as to the true ones, then the errors are given by Eq. (|14p . 

i.e. Ap a = v/t^^ioa' These are represented in the figure 
by the blue dashed lines. 

Figure [16] shows the 2D likelihood surfaces for various 
parameter combinations after marginalizing over all other 
parameters. In all of the panels we take fc max = 1/iMpc -1 
and consider only the 2-cr errors, denoted by A^ 2 = 6.17. 
Again the red solid lines present the results for the case 
where the simulation parameters are fixed and the blue 
dashed lines the case where we marginalize over the simula- 
tion parameters. Clearly, there would be a significant degra- 
dation in the constraining power of any future galaxy clus- 
tering survey, should we not be able to identify the 'true' 
simulation parameters. Note also that there appears to be a 
strong degeneracy between {wo,w a } and {h, fl b }. 
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6.3 Combining information from large-scale 
structures with the CMB 

We now turn to the question of whether adding external 
data sets may help alleviate the degradation of the cosmo- 
logical constraints. Here we only consider the impact on the 
errors of addin g the information from a CMB experiment 
such as Planck (|The Planck Collaborationll2006l '). Note that 
even without Planck data we have already restricted our 
exploration of the cosmological parameter space to include 
only flat models. As described in Appendix fA] we first com- 
pute the CMB Fisher matrix in a set of parameters that are 
suitable for the CMB, and then rotate this matrix to our 
favoured parameter set for describing large-scale structure 
(see also lHilbert et al.ll2012h . We treat the CMB and Large- 
scale structure information as independent and hence the 
Fisher matrices may be added: 



T-Tot _ x-CMB . t-lss 

u fy ft ■ J n ft ~\~ J rv 



*P — ^ap +^o(3 • (32) 

Figure [T7] shows again the errors for the 7 cosmologi- 
cal parameters that we have considered. The differences be- 
tween the errors obtained from marginalizing over the cos- 
mological parameters (red lines) and the cosmological-plus- 
simulation parameters (blue lines) are significantly reduced. 
Thus, inclusion of the CMB information significantly im- 
proves our ability to constrain the cosmological model. 

Figure [TH] shows how the 95%- confidence- level error el- 
lipses changed when we add the CMB information. We see 
again that the constraining power of the combined experi- 
ments significantly improves our ability to constrain cosmol- 
ogy, and also marginalize over the simulation parameters. 



7 CONSTRAINING DARK ENERGY 

We now turn to the question of how well we may constrain 
the time evolution of the dark energy equation of state, i.e. 
{wo, w a }. 

Figure [T5] shows the 95%-confidence-level likelihood 
contours for {wo,w a }, where we use all information from 
LSS up to k = 1/iMpc -1 . The fi gure reveals that the best 



Figure 19. Comparison of the 95%-confidence-level contours in 
the likelihood surface of {wo,w a } after marginalization. The red 
dashed line corresponds to the power spectrum information alone, 
and the green dot-dash line is the same but after marginalizing 
over the simulation parameters. The black solid line is the com- 
bination of the power spectrum information with a Planck-like 
prior, and the blue dotted line is the same but after marginaliz- 
ing over the simulation parameters. 
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Figure 20. Comparison of the dark energy figures of merit. The 
solid red block represents the case of the power spectrum infor- 
mation alone and the green block denotes the same only after 
marginalization over the simulation parameters. The solid blue 
block depicts the case of the power spectrum information com- 
bined with a Planck-like CMB experiment and the magenta block 
denotes the same, after marginalization over the simulation pa- 
rameters. 
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constraints will come from the combination of CMB and LSS 
information. However, if we do not understand how to opti- 
mize our JV-body codes to provide 'true' cosmological power 
spectra, then marginalizing over the simulation parameters 
will be costly for Dark Energy science. 

Currently, the standard way to describe the ability of 
an experiment to cons train w(a) is throug h the figure of 
merit (hereafter FOM) l|Albrecht et al.l l2006). This has been 
defined as the inverse area enclosed by the 2-er error ellipsoid 
of the parameters {wo,w a }. In terms of the Fisher matrix, 
this may be written: 

FOM = 1 , (33) 

7r^/6.17Det[Cov(u;o,™a)] 

where the parameter covariance matrix Cov(ioo,w a ) is the 
2x2 matrix formed from the submatrix of the {wo,w a } 
elements of the inverse of the 7x7 Fisher matrix (e.g. see 
IWang||2008h . 

Figure [20] compares the various dark energy figures of 
merit. From the figure we clearly see that if one marginal- 
izes over the simulation parameters, then there is roughly a 
factor of 2 penalty in our ability constrain the parameters 
{w ,w a }. 



8 CONCLUSIONS 

In this paper we have used a large ensemble of iV-body sim- 
ulations to explore the cosmological information content of 
the matter power spectrum. We have also explored how the 
cosmological information is degraded when we are uncertain 
as to what the 'true' A^-body simulation parameters are. 

In Ej2] we introduced the 'Gemeinsam' likelihood func- 
tion, which takes into account the dependence of the the- 
oretical model on the cosmological and iV-body simulation 
parameters. We reviewed the Fisher matrix formalism for 
forecasting constraints obtainable from measurements of the 
matter power spectrum. The constraints required us to esti- 
mate the fiducial model power spectrum, its covariance ma- 
trix and the first order derivatives of the matter power spec- 
trum with respect to the cosmological and simulation pa- 
rameters. Our fiducial survey consisted of three independent 
volumes, each of which had a volume V M = 3.25 h~ 3 Gpc 3 
but spanning the redshifts z = 0, 0.5 and 1. 

In |3] we described the large ensemble of simulations 
that we have performed in order to compute the Fisher ma- 
trix. We ran 200 simulations to generate the covariance ma- 
trix, we ran 56 simulations to explore the variations of the 
power spectrum with respect to the cosmological parame- 
ters; and 18 simulations to explore its dependence on the 
simulation parameters. 

In |4] we presented the results for the fiducial model. 
We demonstrated that, for k < 0.1/iMpc _1 the errors in 
the power spectrum were reasonably well described by the 
Gaussian prediction. However, on smaller scales the errors 
were found to be significantly larger, and were consistent 
with the presence of a connected trispectrum that scaled as 
Tu oc Pf . We explored the off-diagonal covariance of the 
power spectrum and found that different band powers were 
> 50% correlated for {h,kj} > O.lhMpc -1 at z = 0, and 
for {hi, kj} > 0.2 ftMpc -1 by z = 1. We conclude that non- 



trivial band power correlations must be accounted for in the 
cosmological analysis of large-scale structure data. 

In SJS] we computed the logarithmic derivatives of the 
power spectrum with respect to 7 of the cosmological pa- 
rameters: = {as, fl m , fib, wo, w a , h,n s }. We found that for 
k < 0.1/iMpc -1 the cosmological dependence could be rea- 
sonably well-captured through the variations in the linear 
theory spectra. On smaller scales, the measurements showed 
strong departures from the linear predictions. Interestingly, 
we found that for the parameters {f2 m , fit, n s , h}, and at 
late times, dlog P/da — > 0, as k — > 1/tMpc -1 . This sug- 
gested that there may be very little additional cosmological 
information to be gained on these parameters by the inclu- 
sion of measurements on very small scales. However, we also 
showed that for {a8,wo,w a } the inclusion of small scales 
significantly increases the cosmological information about 
these parameters. 

We then explored the dependence of the matter power 
spectrum on 9 of the simulation parameters used for the 
Gadget-2 code. We found that variations in the choice of 
PMGRID, RCUT, ASMTH, ErrTolForceAcc and Softening could 
in combination lead to percent-level variations in the power 
spectrum on small scales. 

In Sj6] we used the simulations to explore the cosmo- 
logical information content of the matter power spectrum. 
We found that, under the assumptions of flat cosmological 
models, our fiducial survey could constrain {<jg,Q m ,n s } at 
the percent level or better, {fit, h, wo} at the few-percent 
level and w a at the 20% level. However, if we fold into our 
likelihood analysis uncertainties in the simulation parame- 
ters then all of these constraints are degraded by roughly a 
factor of 2. We then showed that adding external data sets, 
such as a Planck-like CMB survey, can help to mitigate the 
effects of marginalization over the simulation parameters. 
In particular, the parameters {n s , h, ttb, 0} are almost un- 
affected by the marginalization procedure. 

In !J7] we focused on the dark energy equation of state 
parameters {wo,w a }- We have shown that marginalizing 
over the simulation parameters significantly degrades our 
ability to constrain the Dark Energy from the power spec- 
trum. Adding the CMB information does help somewhat. 
However, we have computed the dark energy figure of merit 
and found that there is a factor of ~ 2 degradation when 
the simulation parameters are marginalized over. 

In this paper we have worked with the simulation code 
Gadget-2 and a su b-set of paramete rs that are specific to 
it. As discussed in iReed et al.l (|2012h . accurate simulating 
of cosmic structure formation involves more than the code 
parameters. We have neglected to explore the dependence of 
the information on the number of simulation particles, the 
box-size, the initial start redshift. Thus, taken at face value 
this appears to be an optimistic assessment of the problem. 

On the other hand, in principle, a number of the is- 
sues raised in this paper might be mitigated by larger sim- 
ulations: e.g., if one increased the number of particles N 
without limit, then the scale at which the force softening 
modifies the results could be pushed to higher wavenumber, 
since /c so ft 

ee 2Tv/l soCt oc 7V 1/3 . Hence, one could in principle 
find an N sufficiently large that fc so f t will be larger than the 
targeted wavemodes of the designed survey. However, finite 
resources may make this computationally challenging. 

Regarding the generality of our conclusions, one might 
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be concerned that the point in the simulation parameter 
space that we adopted as our fiducial point may bias our 
results and one might ask: how would the results change if 
we adopted another set of fiducial simulation parameters - 
ones closer to the true model? If the likelihood function does 
not vary rapidly over the simulation parameter space, then 
our estimates for the Hessian and and hence the Fisher ma- 
trix will be robust. This, however, is an important question 
which will deserve further attention in the future. We an- 
ticipate that answering it fully will also require one to solve 
the more subtle problem of finding the true set of simulation 
parameters. 

This work has also focused on the problem of simu- 
lating the dark matter only power spectrum. Future work 
will also have to extend this analysis to include the im- 
pact of baryonic physics effects on the clustering due to: 
our approximate handling of the evoluti on of the coupled 
baryo n-CDM fluid after recombination (|Somogvi fc Smith! 
2010); uncertaintie s in the small-scale feedback processes of 
galax y formation |van Daalen et all 1201 ll : ISemboloni et al.l 
1201 ll ). In addition, when exploring alternative cosmologi- 
cal models, new uncertainties will need to be folded into 
the estimates for ex ample if we also wish to constrain the 
dark matter mo del (|Viel et al.l 120121 : ISchneider et al. I l2012l : 
iBird et alJlanl ). 
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The covariance matrices for these observables are: 



12 2 

Cov [Ce,TT, CV.tt] = t — n/J 1 [C(,tt + A^tt] ; (A3) 



./sky It + 1 

Cov [C^TT* C^,ee] = / 



/sky 2£ + 1 



Cf,TE i 



(A4) 



1 2 

COV [CV,TT, CV,Te] = —, 1 ft,TE [C<,TT + A^tt] 



/sky 11 + 1 



COV [CV,EE, CV,Ee] 



/sky It + 1 



(A5) 

[CV.be + Ni,ee\* ; (A6) 



1 2 

COV [CV,EE, C^.te] = 7 oe i C<,TE [C£,EE + A^ee] 



/sky It + 1 



1 1 r 2 

Cov [Cf.TE, CV,te] = — — nn , ; |CV,te 



(A7) 



/sky 2^ + 1 

(CV,EE + A^ > tt)(C , £ i ee + A^ee)] 



(A8) 



APPENDIX A: PLANCK FISHER MATRIX 
Al Computing the CMB matrix 

In computing the Planc k Fisher matrix we f ollow the 
methodology described in lEisenstein et alj (1 19991) and fo r 
the specific implementation we follow lTakada fc Jainl (|2009l V 
We thus assume that the CMB temperature and polarization 
spectra constrain 9 parameters, and for our calculations we 
set their fiducial values to be t hose from the recent WMAP7 
analysis (|Komatsu et al.ll201ll ). The fiducial parameters are: 
dark energy EOS parameters wo = —1.0 and w a = 0.0; 
the density parameter for dark energy Ode = 0.7274; the 
CDM and baryon density parameters scaled by the square 
of the dimensionless Hubble parameter w c = Q c h 2 = 0.1125 
and w b = tt h h 2 = 0.0226 (h = Ho/IlOOkms^Mpc- 1 ]); 
and the primordial spectral index of scalar perturbations 
n s = 0.963; the primordial amplitude of scalar perturba- 
tions As = 2.173 x 10~ 9 ; the running of the spectral index 
a — 0.0; and the optical depth to the last scattering surface 
r = 0.087. Hence we may write our vector of parameters: 



In the above / s k y is the fraction of sky that is surveyed and 
usable for science, and we take / s k y = 0.8. The terms A^tt 
and A^ee denote the beam-noise in the temperature and 
polarization detectors, respectively. These can be expressed 
as: 



A^tt = [wttWfL™^)] 

A^EE = [wEEWi^i^V 1 



(A9) 
(A10) 



where iutt = [A T 9Bcam] 1 and wee = [A B ^Beam] X . The 
beam window function has the form: 

WjcamM = GXp [-£(£ + l)cr B cam] ! CTBeam = 7=P= ■ 

V 8 log 2 
(All) 

For the Planck experiment we assume that we have a single 
frequency band for science (143 GHz), and for this channel 
the f ollowing parameters apply jThe Planck Collaboration! 
2006): angular resolution of the beam #Beam = 7.1' 
[FWHM] ; the beam intensity is A T = 2.2 (T C mb/1K) \pK], 
Ae = 4.2 (Tcmb/IK) [/iK]. We take the temperature of the 
CMB to be T = 2.726K. 



p = (w , W a , f^DE, W c , Wh, T, U B , A s , a) 



(Al) 



The CMB Fish er matrix can be written as 
(|Eisenstein et alJll999h : 



J 7 , 



PaPfj 



EE 



dCe ' x Cov-i[C e , x ,C e , Y ] dCl < Y 



dp,3 



, (A2) 



where {X, Y} G {TT, EE, TE}, where C e ,TT is the tem- 
perature power spectrum, CV,ee is the E-mode polarization 
power spectrum, C^te is the temperature-E-mode polar- 
ization cross-power spectrum. We have been conservative 
and assumed that there will be no significant information 
from the C^bb, the B-mode polarizatio n power spectrum . 
We compute all CMB spectra using CAMB |Lewis et alj|2000h 
and use the additional modu le for time evolving dark energy 
models (|Hu fc Sawicki|[2007] ). We include information from 
all multipoles in the range (2 < t < 1500). 



A2 Transforming from CMB to large-scale 
structure parameters 

In the formation of the large-scale structure we have con- 
sidered how the matter power spectrum depends on the 7 
cosmological parameters: — {erg, fi m , Q n , wq, w a , h, n„} pa- 
rameters. Let us rewrite our original 9-D CMB parameter 
set in terms of a new 9-D large-scale structure parameter 
set. Let us therefore consider the transformation: 



p = (wo,Wa,^E>E,W c ,W h ,T,n B ,A B ,a) 

q = (w ,w a ,Q m ,h,Q b ,T,n s ,a$,a) T . 



(A12) 
(A13) 



Five of the parameters are unchanged from the original set. 
The remaining four are related to the original parameters in 
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the following way: 

Mm = 1 - Ode ; (A14) 
fl b = Mb (1 - 12 DB ) ; (A15) 

»b + We 



(Ts 



^Ap(fc|p)M/(fc7?) 2 



1 1/2 

(A17) 



where P(fc|p) is the matter power spectrum, which depends 
on parameters p a , and where the real space, spherical top- 
hat filter function has the form Wk(y) = 3 [sin y — y cos y]/y 3 , 
with y = kR and R = 8 /i _1 Mpc. 

It can be shown that a Fisher matrix in one set of suit- 
able variables may be represented in another basis space 
through the transformation: 

■^Vto) =^2^c-,J r c,0(p)A 0S (A18) 

where A a7 = dp a /dq~, is the matrix formed from the partial 
derivatives of the old parameters with respect to the new 
ones. From Eqns (|A14|) - (|A17|I we have q a — G a (p), how- 
ever in order to perform the partial derivatives we actually 
require the inverse of these relations, i.e. p a = Ga 1 (q). In 
some cases these inverse relations may easily be determined, 
e.g. Eq. (|A14|I . However, in other cases no analytic inverse 
exists, e.g. Eq. ()A17|) . A simple way around this problem is 
through recalling the following: 

at a 

Hence, if we first compute dq a /dp v = A" 1 ^, then we may 
obtain A, from the fact that: A — > [A -1 ] -1 . 

Let us therefore form the matrix A -1 . For the case of 
those parameters that are unchanged A~^ = S^f}- However, 
for the remaining ones, we have: 

" ilm -1 ; (A20) 



<9f2DE 



dQoE w c + w b 



dw c [w c +wb] 2 
8Q.b w c 



(A21) 

(1 - n DE ) ; (A22) 



dwb [w c + Wb] 



-(I -fi DE ); (A23) 



dh 1 / w c + Wb 



9fi DE 2V(i-n DE ) 



(A24) 



^ = i [(1 - n DE )(w h +w + c)]- 1 ' 2 ■ (A25) 
da a 1 daj _ 1 f d 3 k 8P(k\p) 2 



dp a 2<r 8 dp a 2a 8 J (2ir) 3 dp a 

(A26) 

Note that in order to compute the derivatives dag/dp a we 
use the package CAMB. Besides the generation of various 
CMB power spectra, this package can output the present 
day linear theory matter power spectra P(k\p). The deriva- 
tives are then determined numerically using the standard 
estimator for two sided derivatives. The numerical inverse 
of the mat rix A -1 can easil y be computed using the SVD 
algorithm (|Press et alj|l992h . 
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